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Abstract 

We demonstrate the topological properties of the band-gap of armchair graphene nanoribbons 
in a spatially varying staggered sublattice potential. Several general scaling laws are presented to 
quantify the band gap variation. It is found that all armchair nanoribbons are described by one 
of three distinct classes depending on their width, one of which is the well known massless Dirac 
condition, and the other two we call potentially gapless, and gapless-superlattice. We construct 
an effective theory which faithfully reproduces these results, and makes explicit the nature of the 
competing masses and overlap integrals across a particular sample. Finally we propose several 
systems on which these results should shed considerable light, and which have all already been 
experimentally realized. 

PACS numbers: 73.22.Pr, 73.21.Hb 
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The search for states with topological origins is not a new one. The domain wall states in 
polyacetylene jl] and their field theoretic counterparts j^j], as well as the TKNN invariant and 



the quantum Hall effect j^] followed by the fractional quantum Hall effect and Chern-Simons 
theories are condensed matter examples which each span several decades. Undoubtedly 
this ubiquitous field has recently enjoyed a renaissance due to the excitement surrounding 
topological insulators: time reversal invariant bulk insulators with a single gapless surface 
Dirac cone jl]]. 

The community's interest having been piqued, it is not surprising that graphene, the 
canonical example material with a 2D massless Dirac cone, not to mention the first material 
to have been considered as a topological insulator p, has been cut, bent, layered, crumpled, 
adsorbed, twisted and stretched in order to demonstrate new and exotic topologically non- 
trivial states been a remarkably successful task, and graphene continues to 
produce new and unexpected physics each time someone pokes or prods it in a new way. 

The purpose of the current work is to do this once again. It is well known that quasipar- 
ticles in armchair graphene nanoribbons are, in general, massive. The two K points in bulk 
two dimensional graphene are projected onto each other in momentum space when the sys- 
tem is reduced to a quasi-one dimensional one, allowing the two massless cones to mix, thus 
producing a gap [lsj]. In the seminal paper by Brey and Fertig, this is easily seen by noting 
that the boundary conditions of the two sublattices must match in this special geometry 
[yj]. To contrast, in the case of zig-zag ribbons only one sublattice constitutes each edge, 
thus preserving the gapless K point for all zig-zag edged ribbons. For this single crucial 
reason, zig-zag edged ribbons, samples, and interfaces are the usual candidate systems in 
which topologically borne states are proposed to exist [15] . 

Here, we show that armchair edged ribbons with a staggered sublattice potential have 
a tunable mass-gap which, in the case of gapped ribbons can be decreased, or even closed, 
by varying the potential strength and gradient. In complete analogy with the occurrence of 
massless Dirac ribbons, we show that with increasing width, each effect is observed in every 
third ribbon. By constructing an effective theory which reproduces all tight-binding results, 
we show that this tunability arises from the separate, competing masses originating from 
the staggered sublattice potential and the quasi-one dimensional confinement. This effective 
theory further betrays an analogy with the quantum anomalous Hall effect in pseudo-spin 
space. We propose three realistic systems where these results may be observable: twisted 
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FIG. 1: A typical lattice structure (W/2=18, iy/=7). At the domain wall, the position of the 
alternating positive (red) and negative (blue) sublattice potentials is changed from A+/B- to A- 
/B+, where A and B denote the index of individual pairs. The AB pairs are chosen so as to be 
positioned along the axis of the ribbon. The unit cell is shown by the dashed line. 

bilayer graphene, chemically adsorbed single layer graphene, and lattice or substrate mis- 
matched graphene. 

The Hamiltonian which is the focus of the current work is H = H Q + V, where H is the 
tight-binding Hamiltonian of graphene which is given by 

2W 3 
» 3 

Where ~ 3eV is the overlap integral of first nearest neighbouring sites and is zero 
otherwise, and j denotes the three nearest neighbour vectors. In this formalism, the width 
of a ribbon or superlattice unit cell is W-\/3|5|/2. To create a ribbon we simply construct 
the appropriate unit cell and set t = when a first nearest neighbouring site leaves the edge 
of the ribbon, and to create a superlattice of ribbons with alternating sublattice potential, 
we adopt periodic boundary conditions instead. We also introduce a staggered sublattice 
potential across the width of the ribbon given by 

W/2 

V = s g n (* - l w) m (4°i - C \ +Sl c i+s\ ) ( 2 ) 

i 

Where m is the magnitude of the potential and iw marks the site in the unit cell where 
the potential changes sign, amounting to the location of a domain wall. Although this sharp 
change in sign is a specific choice, we stress that none of our results are qualitatively changed 
by a smooth change in sign of the potential. In Fig. [T] we show a typical ribbon that we are 
considering with the domain wall placed at a typical site iw 
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Upon diagonalising these systems, we first ask the following question: is it possible to 
completely collapse the band gap of an armchair ribbon by applying a periodic potential? 
The answer to this question is yes, but only in one sixth of all possible ribbons, namely those 
where (W — 4)/6 G Z, and only when the iw = W/2. We stress here that these ribbons 
are not Dirac ribbons ((W — 2)/3 G Z). Ribbons where (W — l)/6 G Z almost possess 
this quality, however as they have an odd number of pairs in their unit cells, they cannot 
strictly fulfill the requirements of equally sized domains (W is odd, so iw = W/2 ± 0.5), 
even though the effect is identical. Therefore these ribbons display extremely small minima 
of 6bg compared to their non-Dirac counterparts (W/3 G Z). We label this class of armchair 
ribbons potentially gapless (PG), as a tuned staggered sublattice potential has the ability 
to tune the band-gap to zero, or nearly zero. 

We now move on to the armchair superlattice. The system is now two-dimensional, 
and the lattice translation vectors are orthogonal. In general, the transverse component of 
the Hamiltonian is easily constructed as a two-site unit cell with lattice vectors trivially 
determined by iw- The resulting bandstructure behaviour along the transverse direction is 
cos-like. Due to the oscillating gap along the transverse direction, the band-gap position 
deviates from the T point. In particular, the direct band gap position moves between k x = 0, 
k x = 27r/3 and k x = it, where ir is the dimensionless edge of the Brillouin zone. 

At the gamma point, Dirac and PG superlattices are gapped at m = 0, however PG 
superlattices can close the gap at m = A as in the ribbon case. However the third, so-far 
uncategorised ribbon type here differs. This brings us to a definition for the third class 
of ribbons: those for which W G Z, are always gapless at the T point when m = 0. Thus 
these superlattices are likened to armchair Dirac ribbons, in this sense. This is a particularly 
important class of ribbon, as they show that a very weak modulation across a 2D sample can 
close the gap at the gamma point, thus destroying the valley degeneracy of a 2D sheet. (We 
should point out that this is strictly the low m limit, rather than m = 0, as the latter simply 
corresponds to an infinite graphene sheet with zone-folded (incorrect) energy bands). This 
third class of 'ribbon' we call gapless superlattice (GS) ribbons. The gap at the gamma point 
as a function of m is shown in Fig. [3] by the solid black line. We have thus established that 
all graphene ribbons fall into one of the three classes of Dirac ((W — 2)/3 G Z), potentially 
gapless ((W — l)/3 G Z), and gapless superlattice (W/3 G Z). 

The second important point in the superlattice Brillouin zone is 2n/3, whose behavior 
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FIG. 2: Graphs of band gap trends with increasing m for (top to bottom) GS, PG (with even 
number of A-B pairs), and Dirac ribbons. Note the large dip with all ribbons when iyy = 2, which 
is the smallest possible Dirac domain 

with increasing m is shown in Fig. [3] by the red dashed line. At this point, there is a band 
crossing between the top-most valence band, and the lowest conduction band, leading to 
two Dirac points at m = for Dirac and PG superlattices, directly corresponding to the two 
K points of a 2D sheet. As the mass increases, these bands can mix, and the Dirac points 
disappear. However for some Dirac ribbons, the Dirac points move toward the Brillouin 
zone edge, and can fuse into one, before becoming massive, as can be seen in Fig. [3(c) by 
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the green dotted line. 

Finally, before constructing an effective theory of the ribbons, we first briefly explore 
the relationship between the gradient of the g function of m, ie. dtBc/dm of Dirac 

ribbons, and the location of the domain wall, i\y- A quasi-linear scaling law is observed, 
but ribbons where i\y = 3n + 2 and ribbons where iyy = 3n + 3 have the same gradient 
(for a given n). We find three scaling relationships here, one where the iw is placed such 
that the smaller domain itself makes a Dirac ribbon, one where it makes a GS ribbon, 
and one where it makes a PG ribbon. Respectively, the three scaling laws are given by 
de BG /dm\ m=0 = (1 - 277) - aerify, where, a = 0, 1, 2, and 77 = 1/(W + 1) 

To understand the variation of gap with respect to the mass, we introduce a two-species 
low energy effective model for ribbons with Hamiltonian 

H = H L c\c L + H R c R c R + (H LR c[c R + c.c.) (3) 
Where L and R denote the left and right species, with Hamiltonians given by 

Hi(k) = ka x + (sim + k l ^o y + t' { a z , (4) 

and Hlr = ta^. Here k is along the ribbon axis and so is a good quantum number, 
Siin is the mass from equation [2] where Sj = ±1 for i = L, R. The parameters are justified 
as follows. k l n is the gap obtained from the boundary condition in the original notation 
of Brey and Fertig (with the extra index i denoting left or right species), which is simply 
a quantized k y obtained by taking a finite square well boundary condition and admixing 
the two sublattices. t is the overlap of the wavefunctions in the two sub-systems, which 
is given by t = dyil)\(r y )Hil) R (ry) where H is the original Hamiltonian from equation. 
t\ is an extra mtra-species overlap integral caused by the overlap of each species in the 
other subsystem, such that, for example t' L = dyip\ j (r y )H R ilj L (r y ). This latter term is 
interesting both in that it represents inter-species coupling in another system, but also that 
it produces a mass which is orthogonal to the original mass terms Siirt + k n . 

The mass term s R m + k n causes the decrease in the band gap since s R = —1. This mass 
term changes sign when m = k^, and beyond this value the two species have the same sign 
of their a y mass term. 

The lowest eigenvalues of equation |3] at k = give the gap equation such that 
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FIG. 3: The mass dependence of the band gap at the T point, k x = 2ir/3, and k x = n. The 
behaviour at the T establishes the third universal classification of gapless superlattice 'ribbons', 
which are the only armchair superlattices to become gapless as m — > 0. The other two points define 
the high-symmetry gapless points and band-inversion points, which are lifted for sufficiently large 
m. For each system, iw = 8. 

A = 2min(Aj) where A, are the eigenvalues, and these can be easily solved analytically. 
In particular we show three specific results. Firstly, for = = t' L = t' R = 0, we have 
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A(m,t) = 2{m 2 + t 2 f' 2 (5) 

Which is the monotonic increasing case found in Fig. Ws) for the regular massless arm- 
chair ribbons. Secondly, for = finite but still t' L = t' R = 0, we obtain the similar 
result 



A(m,t,k n ) = 2\k n -{m 2 + t 2 ) 1/2 \ (6) 

Which corresponds to massive armchair ribbons with an equal width in each sub-system. 
The vanishing of the intra-species overlap t' terms in this case has been confirmed numeri- 
cally. Finally, we show the gap equation for t' L = t' R and k^ = k^ which is slightly unphysical, 
but is simpler to look at than the general case while still containing all the essential elements: 



A(m, t, k n , t') = (k 2 n + m 2 + t 2 + t' 2 

-2{k 2 n m 2 + k 2 n t 2 + t'H 2 ) 1 / 2 ) 1 ' 2 (?) 

The effect of the non-zero t' terms, which corresponds to unequal widths of the two 
sub-systems, is to produce a gap at m = ^Jk^ + t 2 where A(m, t, k n ) was gapless. 

So we find that the effective model reproduces all the low energy physics of the larger 
armchair ribbon system. 

These results should prove useful in understanding the physics of several experimentally 
realized materials. The first is adsorbed graphene. Graphene has been successfully hydro- 
genated and fluorinated, which has produced some very interesting results fisL \\\. However 
a convincing physical model for either has thus far proven elusive. The results outlined in 
this work are particularly relevant to adsorbed graphene. It has been shown that randomly 
adsorbed graphene will tend toward uniform adsorption. In particular in the case of semi- 
hydrogenated graphene, it has been shown that hydrogen atoms will tend to bond with 
just one of the two sublattices. However, it is extremely unlikely that randomly adsorbed 
graphene will become completely uniformly adsorbed, and instead should form domains, 
just as those proposed here. Although the domain wall configuration presented here is ideal, 
the variation of the mass-gap as a function of domain wall position should prove extremely 
useful for identifying the behaviour of different adsorbed samples. Further, via controlled 
substrate blocking, it is possible that the adsorption could be tunable in different samples. 
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If possible, a two-gap system could be fabricated, even where the smaller gap could be 
vanishingly small. 

Other systems of particular relevance to these results are twisted bilayer graphene, and 
substrate or lattice-mismatched graphene 18 |. The idea behind these systems is that bilayer 
graphene, or graphene atop a substrate, behaves as a single layer of graphene with a gauge 
field [19I . A twisted bilayer, for example, can be described as a single layer of graphene 
with an oscillating mass term (and gauge field which does not contribute to the gap). As 
the mass term changes sign, it forms a network of domain walls over the sample. It has 
been suggested that these domain walls always lie along an armchair edge [13|. Therefore 
this system forms precisely the systems discussed above, where the magnitude of the mass 
depends on the interlayer spacing. Other forms of lattice-mismatched bilayer graphene, as 
well as substrate mismatched graphene behave analogously. 

In conclusion, we identified two geometric classes of armchair ribbon of period three in the 
width other than the usual Dirac class, namely the potentially gapless class, and the gapless- 
superlattice class. We have shown, by constructing an effeective theory of these materials, 
that due to the competing masses of the sublattice potential and the boundary conditions, 
an ordinarily massive ribbon's gap can be closed, and an ordinarily massless ribbon's gap can 
be opened. It is crucial for future experimental investigations, to understand that ribbons of 
any width can be made gapless, and that bilayer systems in particular can be shown to have 
effective single layer behaviour that necessarily falls into one of these three classes. Finally 
we proposed several systems which have already been realized experimentally to which these 
results are particularly relevant, and we believe should be very useful. 
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